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Abstract. Using the Abstract Cauchy-Kowalewski Theorem we prove that 
the 6-family equation admits, locally in time, a unique analytic solution. More- 
over, if the initial data is real analytic and it belongs to H s with s > 3/2, and 
the momentum density uq — uq xx does not change sign, we prove that the solu- 
tion stays analytic globally in time, for b > 1. Using pseudospectral numerical 
methods, we study, also, the singularity formation for the 6-family equations 
with the singularity tracking method. This method allows us to follow the 
process of the singularity formation in the complex plane as the singularity 
approaches the real axis, estimating the rate of decay of the Fourier spectrum. 
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1. Introduction 

(N ; In 2003 Holm and Staley pH [H] introduced, as a one-dimensional version of 

. active fluid transport, the b-family equation 

(1.1) m t + um x + bmu x =0, t > 0, x G E, 

OV 

: 

■ where b is a real parameter, u — u(x,t) the fluid velocity, and m = u — u xx the 
momentum density. The quadratic terms in (|l.ip represent the competition and 

O^l ■ balance, in fluid convection between nonlinear transport and amplification due to 

6-dimensional stretching. The parameter b characterizes the equations because it 
is the ratio of stretching to convective transport and also the number of covariant 
dimensions associated with the momentum density m. 
. Within this family of evolutionary equations we have both the Camassa-Holm 

equation [H 115] in correspondence of b— 2 and the Degasperis-Procesi equation 
[5] in correspondence of b — 3 as special cases. Using the method of asymptotic 
intcgrability, it has been shown that these arc the unique two cases in which (jl.lj) 
is completely integrable [8j. A quite complete list of references on the Camassa- 
Holm, the Degasperis-Procesi equations and the b-family equations can be found 
in H0 [251 111 [23 [28]. 

We can rewrite (jl.l[) is the following third order form in terms only of the 
unknown u 

(1.2) u t - u txx + (b + l)uu x = bu x u xx + uu xxx , t > 0, x G K. 

Moreover, (jl.l[) and (|1.2p are also equivalent to the following hyperbolic elliptic 
system 

b 3 b 

(1.3) u t + uu x + P x = 0, -P xx + P = -u 2 + u\. 

i 
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Since e x /2 is the Green's function of the Hclmholtz operator 1 — d 2 we can rewrite 
(|1.3p in the following integro-differential form 

r e -\x-y\ / b 3 — 6 \ 
u t +uu x + J sign(y-x) i^u(y,t)-\ — u x (y,t) \ dy = 0. 

In this paper we are concerned with the initial value problem associated to 
(|1.2[) . therefore we augment the equation with the initial condition 

(1.4) u(x,0) = u (x), iel. 

Assuming that 

3 

u o ~ UQ yXX has constant sign 

using the Kato's theory and giving a precise description of the blow-up scenario, 
Gui, Liu, and Tian proved the existence of global in time solutions. 

The rest of the paper is organized as follows. In Section |2~T1 we introduce the 
functional spaces in which our results hold. Section 12.21 is devoted to the local 
in time existence of solutions for the Cauchy problem associated to (|1.2p . The 
existence of global in time analytical solution for (|1.2p is proved in Section [3] In 
Section|¥]we analyze the singularity formation using the singularity tracking method 
for different values of b > — 1. 

2. Local in time result 

2.1. Functional settings and the ACK Theorem. In this section we intro- 
duce the functional spaces where local in time analyticity for the 6-family equations 
will be proved (see 17J[18j|22j). We consider x a complex variable, with $lx and 
$sx respectively its real and imaginary part. We consider f(x) a complex function 
analytic in 

(2.1) D{p) ={i£C : \^x\ < p} 

and / G L 2 (T(^sx)) for ^sx G {—p,p). The path where the L 2 integration is per- 
formed is: 

(2.2) T(b) = {x G C : \<3x\ = b} . 

This means that if Qx is inside (—p,p), then f{tftx + i^sx) is a square integrablc 
function of 5Rx. We denote by H°' p the set of complex functions f(x) with the 
above properties and such that 

(2.3) \f\ P = sup ||/(- + i9x)||i y a(r(ste)) < oo 
is satisfied. 

We denote by H k,p the space of all complex functions f(x) such that 

(2.4) <>:/• // :, / '. for < j < k; and ||f|| k>p = £ |^f| p . 

0<j<k 

In term of the Fourier transform /(£) of /, the previous norm can be rewritten as 

(2.5) ii/iu,,= (/ e 2 ^(i + \e) k \m)\ 2 d^j 1/2 ■ 

We denote by Hl% the space of all functions f(x,t) G C 1 ([0,T\,H 1 > p ) and 

\f\t,p,P,T= Sup \f(',t)\ 0tP -pt+ Sup \dtf(;t)\ 0tP -pt+ Sup \d x f(;t)\o, p -pt 
0<t<T 0<t<T 0<t<T 

is finite. In an analogous way, Hg'j, is the space of all functions 
/ G C k ([0,T],H k <P) and \f\ k , P ,p, T < oo. 
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We conclude this section recalling the ACK theorem in the form given by 
Safonov [21]. 

Consider the problem 

u = F(u, t), 

where u 6 X p and X p is a Banach space with norm denoted by | • \ p . 

Theorem 1. ACK Theorem: Suppose that 3R > 0, p Q > 0, and /3 > such 
that if < t < po/ Pq, the following properties hold: 

(1) for all < p' < p < po and u G X p such that sup 0<t<T |u(-,i)| p < R, the 
map F(u,t) : [0,T] H- X p i is continuous; 

(2) for all < p < po the function 

F(0,t):[0,p /f3 }^{ueX p : sup \u(-,t)\ p <R} 

0<t<T 

is continous and \F(Q,t)\ p < Rq < R: 

(3) for all < p' < p(s) < pq and for all w and u € {u € X p : 
su Po<t<T K-^)Ip-/3 « < R }> results 

\F(u,t) - F(w,t)\ p , < C f l "~ W|p( ; ) ^; 

Jo Pis) - p' 

then there exists j3 > /?o such that for all < p < po equation () has a unique 
solution u(-,t) e X p with t £ [0, (po — p)/Po] and moreover 

sn Vp< P0 -Pt H;t)\ P < R- 

2.2. First main theorem: short time existence of analytic solution. 

In this section we proof the following first main theorem: 

Theorem 2. Local in time analyticity: Let the initial data for the b- family 
uq 6 H > p ° . Then there exists ft > such that for any p, with < p < pa, there 
exists a unique continuously differentiable w.r.t. time solution u of the b-family 
equation with the following property: 



• u(; t) £ H l >P and d t u(-,t) € H 1 ' p , when t £ 



n po-p 



Proof. We apply the abstract Cauchy-Kovalevskaya (ACK) theorem to the 
6- family equation to prove the short time existence of analytic solution. We remand 
to [18j for the details and here we give only a sketch of the proof. 

To apply the ACK theorem to the 6-family equation, we write the 6-family 
equation as: 

(2.6) u t + uu x = -A- 2 Qw 2 + ^y^-u 2 x 

where A 2 = (1 — d 2 ), and, with an integration in time, we transform (|2.6I) in 

(2.7) u = F(u,t), with 

/ h i _ h \ 1 

dt! . 



F(u,t) = it 



A' 



b 2 3—6 2 
2 U + — U * 



It is clear that properties 1 and 2 of the ACK theorem are satisfied by the operator 
F(u,t) given in (|2.7I) . For the validity of hypothesis 3, one considers first the 
norm \F(u,t)\i_ p i , where the bound is easily obtained using the following estimates 
which are a direct consequence of the Cauchy estimate on the derivative of analytic 
function (22 [22]): 

Lemma 1. Let u,w e H 1,p " and p' < p" , then \\uw x — wu x \\i tP > < c Jflp/' ■ 

Finally, the bound on \d x F(u, t)\i tP i is obtained using the fact that the operator 
A~ 2 is bounded. □ 
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3. Global in time analyticity 

In this section we prove that the 6-family equation admits a strong solution in 
Gevrey-class functional space of index 1. This analysis gives an explicit lower bound 
on the shrinking of the radius of analyticity of the solution of the ^-family equation 
with given analytical initial data, and assuming suitable regularity in Sobolev space. 

We consider the space D(A r e pA ), which consists of the real restriction of com- 
plex analytic functions on the strip D(p). The operator A r e pA is the Gevrey op- 
erator with A 2 = (1—8%.)- The space D(A r e p ) is equivalent to the space H r ' p 
defined in section |2~T1 with norm given by (|2.5[) . then we define the Gevrey class 
G 1 of index 1 by (see [16]): 

(3.1) G 1 = (J H r ' p . 

p>0 

We suppose that the initial data uq € D(A r e Po ), with po and r positive. 
In this section we prove the following second main theorem: 

Theorem 3. Global in time analyticity: Let b > 1 and A 2 = (I — d 2 ). Let 
uq 6 D(A r e PaA ), with r > 3/2, po > and mo — uo — uq xx does not change sign. 
Then the unique solution u of the b-family equation lies in Gevrey class of index 1 
globally in time. 

Proof. We start the energy type estimate on D(A r e pA ) for the 6-family equa- 
tion. Supposing that p is a decreasing, positive C 1 function of the time t with 
p(0) = po, we consider the L 2 -scalar product of the b- family equation (|2.6I) with 
A 2r - 2 e 2pA u: 

l±\\A r e pA u\\ 2 -p\\A^ 2 e pA u\\ 2 = 

= \(A r e pA uu x ,A r e pA u) + (A r ~ 1 e pA {^u 2 + ^u*) ,A r ~ 1 e pA u)\ 
< \{A r e pA uu x ,A r e pA u)\ + M|(A^V A (u 2 ) x , A r ~ l e pA u)\ 

(3.2) +\^l\ { A^e pA (u 2 x ) x ,A^e pA u)\. 

Now, we start to estimate the three terms on the second member of (|3.2p . 
h = \{A r e pA uu x ,A r e pA u)\, I 2 = \{A r ~ l e pA (u 2 ) x , A r ~ 1 e pA u)\ and 
h = \{A^e pA {ul) x ,A^e pA u)\. 
From Lemma below, we have 

(3.3) Ji < dWA^W 3 + C 2 p\\A r e pA u\\\\A r+1 / 2 e pA u\\ 2 . 
Using Cauchy estimate and Lemma [2] below we obtain: 

(3.4) I 2 < (A r e pA u 2 ,A r - 1 e pA u)\ < \\A r e pA u 2 \\\\A r ~ 1 e pA u\\ 

< \\A r e pA u\\ 2 \\A r - 1 e pA u\\ < \\A r e pA u\\ 3 

< dWA'-uf + pC 2 \\A r+1 ^e pA u\\ 3 

< Ci||A'u|| 3 + P C 2 (\\A r e pA u\\ 1/3 \\A r+1/2 e pA u\\ 2/3 y 

< dWA^f + pC 2 \\A r e pA u\\\\A r+1 ^ 2 e pA u\\ 2 - 
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(3.5) h < (A r - 1 e pA u 2 x ,A r eP A u)\ < \\A r - 1 e pA ul\\\\A r e pA u\\ 

< ll/T-V-VlHlyTe^ull < \\A r e pA u\\ 3 

< dWA'-uf + P C 2 \\A r+1 ^e pA u\\ 3 

< C4A r uf + pC 2 (||A'' e '' A u || 1 / 3 ||^+ 1 /V A U || 2 / 3 ) 3 

< C^u^ + pC 2 \\A r e pA u\\\\A r+1 l 2 e pA u\\ 2 . 
Collecting all the previous estimates (|3.3[) . p.5[) and p.6[) . we obtain 

(3.6) \j t \\ A ^ pAu W 2 < (p + C 2P \\A r e pA u\\) WA^e^uf + CtWA^f, 

where constants C\ and C 2 depend only on b and p. The last inequality (I3.6[) 
implies that u g D(A r e pA ) with 

(3.7) p(t) = poe - c ^o^ ArePoA ^+ c ^fo II^H 8 *"))*'. 

To prove the global in time analyticity of the solution u, it is necessary to have 
a priori Sobolev estimate on ||^4 r u||. To do that we use the global existence result 
for the 6-family equation, given by the result of Lemma ([][]) below. □ 

Lemma 2. Let p > and let u,v e D(A r e pA ) with r > 1/2; and let w e 
D(A r+s e pA ), with p,r, s > 0. There exists constants C, C\ and C 2 , depending only 
on r, such that 

\\A r e pA (uv)\\ < C\\A r e pA u\\\\A r e pA v\\; 
\\A r e pA w\\ < CiP r HI + pC 2 \\A r+s e pA w\\. 
where \\ ■ | is the usual L 2 norm. 

Proof. See [16]. □ 

Lemma 3. Let v e D(A r e pA ) and u G D{A r+1 / 2 e pA ), with p > and r > 3/2, 
then there exists constants C\ and C 2 , depending only on r, such that 

\(A r e pA vu x ,A r e pA u)\ < d \\A' 'v\\ \\A r u\\ 2 + C 2 p\\ A r e pA v\\ \\A r+1 / 2 e pA u\\ 2 . 

Proof. See [1] and [18j. □ 

Lemma 4. Let b > 1 and let u E H S (R) n L lr R), with s > 3/2 and mo = 
uo — uqxx does not change sign, then the corresponding solution does exists globally. 

Proof. See pj] for the proof. □ 

4. Singularity formation for the 6-family equation 

In this section we present the analysis of the singularity formation for the b- 
family equation using the singularity tracking method for different values of the 
parameter b. The singularity tracking method is based on the following Laplace 
asymptotic formula (see [5]) for the Fourier spectrum of a complex function u(z): 

(4.1) u k - C|/cr (1+a) exp(-<5|fc|) exp(i\k\x*), 

where z* = x* + iS is the nearest complex singularity to the real axis, a is the 
algebraic type, and at the singularity u(z) ~ (z — z*) a . Therefore an estimation 
of 5 and x* gives the complex location of the singularity. Resolving the rate of 
algebraic decay 1 + a, one can classify the singularity type. In particular, when 
one deals with a PDE with solution u(t), if at a critic time t s results that 5(t s ) =0 
then the width of analyticity of the solution is zero and one can observe a blow 
up of the solution or a blow up of its derivatives. The reader is invited to see 
[II 11 HH [101 [201 [M [21] for more details on this method. 
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To apply the singularity tracking method to the 6-family, we numerically solve 
the equations in the spatial domain [— tt,tt], employing periodic boundary condi- 
tions. We write the 6-family equations in pseudodiffcrcntial form 



Ui 



= -A~ 



3-6 



2 2 

and the dynamic of the kth Fourier mode of u is described by the following ODE 

ik ( b r^Ts 3 — b - 



(4.2) 



d t u k = - 



(uu x ) 



k 



1 



b r-^r 3 

2 (U) ^ + — 



'-(ul) 



Dividing the time interval [0,T] in N sub-intervals of size At = T/N, we approxi- 
mate the solution u as 

K/2 
k=~K/2 



u(x, nAt) 



Akx 



and we solve the system of ODEs (|4.2p using explicit Runge-Kutta method of the 
4th order with initial conditions given by 

u° = u 0k , 

where uo k are the Fourier coefficients of the initial data u(x, 0) = uq(x). The reader 
is referred to [9] for the analysis of the convergence properties of this numerical 
scheme for the Camassa-Holm equation case, which is the ^-family equation with 
b = 2. 

The estimation of the parameters C, a and 5 in the Laplace asymptotic formula 
(14. ip is performed with the technique of the sliding-fitting with length 3 (see [3l 
[9l 1191 1201 123| for details) . This procedure is based on the local estimation of the 
wanted values for each k mode, by the formulas: 



lO! 



(4.3) 

(4.4) 
(4.5) 



a{k) 

8{k) 
logC(fc) 



V K 



log 



log 



k 2 

(fe-l)(fc+l) 



Uk 



a log 



u k+ i J \k + l 

log |u fc | + a(k) log (k) + kS(k). 



The asymptotic values for large k can be computed with an extrapolation process 
using the epsilon algorithm of Wynn (see |19]V giving the wanted values C, a and 
6 . 

For all the numerical simulations we shall present in this work, we have used 
up to 4096 Fourier modes calculated with 32-digit (using the ARPREC package [2]). 

4.1. Numerical solutions. For the numerical results we consider two differ- 
ent initial data: the first one is uq — sin(x) and the second one is uq = 1 + sin(:r), 
and in the sequel we call them type I e II respectively. These initial conditions does 
not satisfy the hypothesis of theorem 2, therefore the solutions arc candidate to 
develop a singularity in a finite time. This work want to extend the results found 
in [9] , where the singularities formation for Camassa-Holm equation was analyzed. 
In particular, for initial condition of type I, Camassa-Holm equation develops a 
singularity in a finite time of algebraic type 3/5, weaker with respect to the cubic- 
root singularity developed by Burgers equation |24j . probably due to the presence 
of the dispersive terms. For initial condition of type II, Camassa-Holm equation 
develops in a finite time singularity of algebraic type 2/3, which is the behavior of 
the peakons solution for Camassa-Holm equation. 
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Here we present the numerical results of the Degasperi-Procesi equation, which 
corresponds to the 6-family equation with b — 3 with respect to initial data of type 
I and II. 



b=3 u Q =sin(x) 




Figure 1. The behavior in time of the numerical solution of 
Degasperi-Procesi equation with initial condition of type I. At time 
t s = 0.8295 a singularity forms as a blow up of the first derivative. 




Figure 2. The behavior in time of the numerical solution of 
Degasperi-Procesi equation with initial condition of type II. At 
time t s = 0.8875 a singularity forms as a blow up of the first de- 
rivative. 

The behaviors in time of the numerical solutions and their derivatives of the 
Degasperi-Procesi equation with initial condition I and II, are shown respectively 
in FigfT] and FigJ5] The critical times when singularities develop are respectively 
t s = 0.8295 and t s — 0.8875. At these critical times, it is evident a blow up of the 
first derivatives of the solutions. 

The analysis of the singularities formation and their characterization is re- 
manded in the sequel, where we perform a complete analysis for all the 6-family 
equations. 

5. Classifications of complex singularities: initial data uq = sin(a;). 

In this section we report the results obtained by applying the singularity track- 
ing method to the 6— family equation with initial datum of type I. 

In Fig |3] it is shown the behavior in time of the spectrum of the Degasperi- 
Procesi numerical solution (6 = 3), starting at time t = 0.35 up to singularity time 
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t s = 0.8295 with increments of 0.05: at critical time t s the spectrum totally loss 
exponential decay, as a singularity of cubic-root type hits the real axis, with the 
blow-up of the first derivative of the solution. The extrapolation of S(t) and a(t) 
(on Fig|4]) is obtained by the sliding fitting procedure of length 3 




O 500 1 OOO 1 500 2000 



Figure 3. The behavior of the spectrum in time of the Degasperi- 
Procesi equation numerical solution with initial datum of type I, 
starting at time t = 0.35 up to singularity time t s — 0.8295 with 
increments of 0.05. 



t =0.82 b=3 u =sin(x) 




Figure 4. The tracking singularity method results of the 
Degasperi-Procesi equation numerical solution with initial datum 
of type I. On the top, the evolution in time of the width of the 
analyticity strip 5. On the bottom, the evolution in time of the 
algebraic character a. The algebraic character at singularity time 
is a(t s ) = 0.33. 

We have also performed several other numerical simulations for the initial con- 
dition I, by varying the parameter b > — 1. We have obtained that in each case a 
singularity develops in a finite time, and in FigJS]and FigJB]the algebraic character 
of the singularities and the times of singularity formation are shown with respect to 
the parameter b. We remember that for b = — 1 the solution for the initial datum of 
type I is a traveling wave solution, and therefore in this case the solution is globally 
analytic in time. 
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u =sin(x) 



Figure 5. The times of singularity for the b- family equation with 
initial condition of type I with respect to the parameter b. On the 
right, the singularity algebraic characters a for the b- family equa- 
tion with initial condition of type I with respect to the parameter 
b. 



2.5 r 

«(t ) 



u n =sin(x) 



-1 o 



Figure 6. The singularity algebraic characters a for the b- family 
equation with initial condition of type I with respect to the param- 
eter b. 



6. Classifications of complex singularities: initial data uq = 1 + sin(x). 

In this section we report the results obtained by applying the singularity track- 
ing method to the 6-family equation with initial datum of type II. 

In FigJ7]it is shown the behavior in time of the spectrum of of the Degasperi- 
Procesi numerical solution, starting at time t = 0.4 up to singularity time t s — 
0.8875 with increments of 0.05. On Fig|S]the evolutions in time of the width of the 
analyticity strip S and the algebraic character a are shown. One can see that at 
the singularity time t s — 0.8875 the algebraic character is a{t s ) = 0.41, meaning 
that a blow-up of the first derivative of the solution occurs. The results are also 
obtained by applying the sliding fitting procedure with length 3. 

Also for the initial condition of type II we have performed several other numer- 
ical simulations by varying the parameter b. The singularity time and the algebraic 
character of the singularities are shown in Figj9]and Fig|10]w.r.t the parameter b. 
As observed in the case of the initial condition of type I, for b = — 1 the solution for 
the initial data of type II is a traveling wave solution, and therefore the solution is 
globally analytic in time. 
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k 



O 500 1000 1500 2000 



Figure 7. The behavior in time of the spectrum of the Degasperi- 
Procesi equation numerical solution with initial datum of type 
II, starting at time t = 0.4 up to singularity time t s — 0.8875 with 
increments of 0.05. 





Figure 8. The tracking singularity method results of the 
Degasperi-Procesi equation numerical solution with initial datum 
of type II up to singularity time t s — 0.8875. On the top, the 
evolution in time of the width of the analyticity strip S. On the 
bottom, the evolution in time of the algebraic character a. The 
algebraic character at singularity time is a(t s ) — 0.41. 



7. Conclusions 

In this paper we have considered the ^-family equations, introduced by Holm 
and Staley [131 114] , Our main interest has been to show the well-posedness of 
the 6-family equation in analytic function spaces. Using the Abstract Cauchy- 
Kowalewski Theorem, we prove that the 6-family equation admits a unique analytic 
solution, both locally in time (Theorem (|2|) then globally in time (Theorem ([3])). 
An interesting problem in the theory of PDEs is the singularity formation and how 
to detect it from a numerical point of view. One of the most effective method 
to achieve this goal is the singularity tracking method. It consists to follow the 
singularity in the complex plane before the appearance as a blow up of the solution. 
The distance of the singularity from the real axis can be estimated through the study 
of the exponential rate of decay of the Fourier spectrum. The algebraic character 
of the singularity is determinate, instead, by the the study of the algebraic rate 
of decay of the Fourier spectrum. We use numerical pscudospectral methods to 
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u =1 -i-sin(x) 



Figure 9. The times of singularity for the b- family equation with 
initial condition of type II with respect to the parameter b. 



oc(t ) 



u =1 -t-sin(x) 



Figure 10. The singularity algebraic characters a for the b- family 
equation with initial condition of type II with respect to the pa- 
rameter b. 

study the process of singularity formation for the 6-family equations, w.r.t the 
parameter b > — 1, and the singularity time and the algebraic characters are shown 
in Sections [5] and \5\ We analysed also the cases for 6 < —1 and also in these 
cases the solutions develop a singularity, but it was complicated to determinate 
their algebraic character. Very high numeric precision is needed and it will be the 
subject of future work. 
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